METHOD FOR EVALUATING COMPLEX REFRACTIVE INDICIES UTILIZING IR 



RANGE ELLIPSOMETRY 

This Application is a CIP of Provisional Application Serial No, 
5 60/212,848 filed 06/26/00, and of Utility Application Serial No. 
09/474,318 filed 12/29/99. 

TECHNICAL FIELD 

10 The present invention relates to characterisation of 

optically thick material systems, and more particularly is a 
method of simultaneously evaluating mathematical model parameters 
,^ including Euler angles, and directions and magnitudes of complex 
|j components of orthogonally related dielectric functions or 
l|j refractive indicies in an optically thick material system which 
g presents with an optical axis oriented either in-plane or 
y| out-of-plane with respect to an alignment surface thereof. The 
present invention method, while not limited thereto, is 
particularly applicable in investigating optically thick material 
2(fj systems which are at least uniaxial in that magnitudes of 
HI corresponding real, and magnitudes of corresponding imaginary, 
yj components of at least two orthogonally related optically thick 
71 material system characterizing diagonalized tensor: 
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complex dielectric functions or refractive indicies are equal. 
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BACKGROUND 



Ellipsometry is a well known means by which to monitor 
optical and physical properties of material systems, (eg. a 
substrate which possibly comprises thin films on a surface 
thereof), in brief, ellipsometry provides that a polarized beam 
of electromagnetic radiation of one or more wavelengths is caused 
to impinge upon a material system along one or more angles of 
incidence and interact, (eg. reflect from or transmit through), 
therewith. Beams of electromagnetic radiation can be considered 
as comprised of two orthogonal components, (ie. "P" and "S"), 
where "P" identifies a plane which contains both an incident beam 
of electromagnetic radiation, and a normal to an investigated 
alignment surface of a material system being investigated, and 
where "S M identifies a plane perpendicular to the "P" plane and 
parallel to said alignment surface of said material system. A 
change in polarization state in a polarized beam of 
electromagnetic radiation caused by said interaction with a 
material system, is representative of properties of said material 
system. (Note that Polarization State basically refers to a 
magnitude of a ratio of orthogonal component magnitudes in a 
polarized beam of electromagnetic radiation, and a phase angle 
therebetween, although absolute values for orthogonal components 
and "handedness", and even percent of polarization are further 
full polarization state determining factors). Generally two well 
known angles, (PSI and DELTA), which characterize a material 
system at a given Wavelength and Angle-of -Incidence, are 
determined by analysis of data which represents change in 
polarization state. PSI is a ratio of the "P" and "S" component 
magnitudes, and DELTA is the phase angle therebetween. 

For general interest, it is noted that spectrophotometer 
systems do not provide data regarding DELTA related phase angle, 
and often provide absolute, rather than PSI {'jl ) type ratio, 



intensity values, For instance, utilizing a transmission Jones 
matrix: 
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to represent mathematically how a material system affects a 
polarized electromagnetic beam, with Epi and Esi orthogonal 
components, and which is caused to interact therewith, it is to 
10 be appreciated that an ell ipsometer/polar imeter system might 
return an on-diagonal ratios 

(Tpp/Tss) 

lyiS and off-diagonal ratios 

=2 (Tsp/Tss) 

A (Tps/Tss ) = 

201 

(Tsp/Tpp) = 

% (Tps/Tpp) = 

25 as a function of a selection from the group consisting of: 

angle -of -incidence; 
wavelength; and 

material system rotation about a normal to 
2q said alignment surface ♦ 

A spectrophotometer system might return Tpp, Tss, Tsp and/or Tps 
as a function of angle-of-incidence; wavelength; and/or material 
system rotation about a normal to said alignment surface. It is 
to be understood that a reflection Jones Matrix 
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Epoj = jRpp Rspj 
JEso| = [_Rps Rssj 

could have also been utilized for said demonstrat ion, and that 
5 generally the present invention method can be practiced utilizing 
data obtained with an ellipsometer system configured in 
transmission or reflection modes, and/or combinations thereof in 
any steps thereof, wherein functionality is preserved. 

10 Continuing, Ellipsometer Systems generally include a source 

of a beam of electromagnetic radiation, a Polarizer, which serves 
to impose a linear state of polarization on a beam of 
electromagnetic radiation, a Stage for supporting a material 
p system, and an Analyzer which serves to select a polarization 
ISO state in a beam of electromagnetic radiation after it has 
m interacted with a material system, and prior to being passed to a 
CO Detector System for analysis. As well, one or more 
^ Compensator ( s ) can be present and serve to affect a phase angle 
38 between orthogonal components of a polarized beam of 
20L electromagnetic radiation- This is especially important where 
§1 it is necessary to determine the "Handedness" of a polarized beam 
Hi 0 f electromagnetic radiation . 

ftp 

i'"~ r i 

A number of types of ellipsometer systems exist, such as 
25 those which include rotating elements and those which include 

modulation elements. Those including rotating elements include 
Rotating Polarizer (RP), Rotating Analyzer (RA), Rotating 
Compensator (RC) . It is noted that Rotating Compensator 
Ellipsometer Systems do not demonstrate "Dead-Spots" where 
30 obtaining data is difficult. They can read PSI and DELTA 

(A) of a material system over a full Range of Degrees with the 
only limitation being that if PSI becomes essentially zero (0.0), 
one can not then determine DELTA as there is not a sufficient PSI 
([jJ ) Polar Vector Length to form the angle between the PSI ( \p) 



Epi 
Esi 



4 



Vector and an "X" axis. In comparison, Rotating Analyzer and 
Rotating Polarizer Ellipsometers have "Dead Spots" at DELTA'S (A) 

near 0.0 or 180 Degrees and Modulation Element Ellipsometers 
also have "Dead Spots" at PSI ({// ) near 45 Degrees. The present 
invention method can be practiced with essentially any 
ellipsometer system. 

A Search of Patents relevant to the present invention has 
identified very little of specific relevance. One Patent, to 
Dill, No. 4,053,232 describes a Rotating-Compensator Ellipsometer 
System, which operates utilizes monochromatic light. Two Patents 
which identify systems which utilize Polychromatic light in 
investigation of material systems are described in Patent Nos. 
5,596,406 and 4,668,086, to Rosencwaig et al. and Redner, 
respectively, were also identified. Also identified is a Patent 
to Woollam et al, No. 5,373,359 as it describes a Rotating 
Analyzer Ellipsometer System which utilizes white light. Patents 
continued from the 359 Woollam et al . Patent are, Nos. 5,504,582 
to Johs et al. and 5,521,706 to Green et al. Said 582 Johs et 
al. and 706 Green et al. Patents describe use of polychromatic 
light in a Rotating Analyzer Ellipsometer System. A Patent to 
Bernoux et al., No. 5,329,357 is identified as it describes the 
use of optical fibers as input and output means in an 
ellipsometer system. A Patent to Chen et al., No. 5,581,350 is 
identified as it describes the application of regression in 
calibration of ellipsometer systems. 

A particularly interesting Patent to Herzinger, No. 
5,835,222, is identified, and incorporated hereinto by reference, 
as it describes an ell ipsometr ic based method for identifying the 
orientation of an optical axis in a material system with respect 
to an alignment surface thereof. One recitation of the method of 
the 222 Patent method of qualifying a material system as having 
an optical axis oriented in a desired locus with respect to an 



alignment surface thereof, comprises, in a functional order, the 
steps of: 

a. by ellipsometric techniques determining the magnltude(s) of 
at least one member of the group consisting of: 

real; imaginary; and a combination of real and 
imaginary components; 

for at least one ratio of components of a material system 
representing transmission Jones matrix as a function of at least 
one member of the group consisting of: 

wavelength; and 

"P" plane angle-of -incidence of an investigating 
polarized beam of electromagnetic radiation to 
an alignment surface of said material system; and 

optical axis radial direction rotation angle; 

said at least one ratio of material system representing 
transmission Jones matrix components being selected from the 
group consisting of: 

on-diagonal ratio 

(Tpp/Tss } ; 

and off-diagonal ratios 

(Tsp/Tss) ; 
(Tps/Tss ) ; 
(Tsp/Tpp) ; 
(Tps/Tpp) ; 

where Tpp, Tss, Tps and Tsp are the components of a two 
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dimensional material system representing transmission Jones 
matrix: 



Epoi = 


Tpp 


Tsp 


Epi 


i ! 
: Eso] = 


Tps 


Tss 


Esi 

_ — > 



which material system representing transmission Jones matrix 
describes the change in polarization state between said incident 
and transmitted portions of said beam of electromagnetic 
radiation; 

b. providing a mathematical model of said material system 
comprising at least one deviation angle calibration parameter 
which represents a non-coincidence of said optical axis with a 
desired locus with respect to said alignment surface, said 
mathematical model serving to relate indices of refraction, 
thickness and optical axis direction over a range of at least one 
member of the group consisting of: 

wavelength; and 
f, P" plane angle -of -incidence of an investigating 
polarised beam of electromagnetic radiation to 
an alignment surface of said material system; and 
optical axis radial direction rotation angle; 

c. performing a regression procedure of said mathematical 
model onto said magnitude of at least one member of the group 
consisting of: 

real; 

imaginary; and 

a combination of real and imaginary components; 



of at least one step a. selected ratio of material system 



representing transmission Jones matrix components as a function 
of at least one member of the group consisting of; 

wavelength; 

"P" plane angle-of -incidence of an investigating 
polarized beam of electromagnetic radiation to 
an alignment surface of said material system; and 
optical axis radial direction rotation angle; 

such that said at least one deviation angle calibration parameter 
is evaluated * 

Material systems are accepted if the deviation angle is small 
enough to fall within selected guidelines. 

Said method of qualifying a material system as having an optical 
axis oriented in a desired locus with respect to an alignment 
surface thereof can further comprising the step of; 

e* plotting magnitude(s) of at least one determined on-diagonal 
(Tpp/Tss) and/or of f -diagonal ratio(s) selected from the group 
consisting of: 

(Tps/Tpp) and (Tsp/Tss) and (Tps/Tss) and (Tsp/Tpp); 

with respect to at least one parameter selected from the group 
consisting of: 

wavelength; 

"P" plane angle-of -incidence of an investigating 
polarized beam of electromagnetic radiation to 
an alignment surface of said material system; and 
optical axis radial direction rotation angle; 

and accepting said material system based upon observing said 
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plot(s) only if said plot(s) essentially match known "template" 
plots . 

Another Patent, No. 5,757,494 to Green et al., is identified 
as it describes a method for enabling investigation of a material 
system with ellipsometer systems, even in ellipsometr ic PSI 
and/or DELTA regions wherein said ellipsometer system is, without 
the invention, "blind", (eg. DELTA of (0.0) or (180) degrees) in 
Rotating Polarizer and Rotating Analyzer Ellipsometer Systems, 
and/or (PSI or (45) degrees in Modulation Element Ellipsometer 
Systems). Said 494 Patent describes a method of determination of 
material system PSI and DELTA values with improved accuracy and 
precision comprising utilising a variable compensator, said 
method comprising in a functional order, the steps of: 

a. providing an ellipsometer system selected from the group 
consisting of: 

modulation element, rotatable element and rotating element; 

which ellipsometer system enables accurate and precise 
determination of PSI and DELTA values of essentially any 
investigatable material system; said ellipsometer system 
comprising means for setting at least one polarization state in a 
beam of polarized light and means for identifying, and means for 
monitoring a polarization state in said polarized beam of light, 
after an interaction thereof with a material system; 

between said means for setting at least one polarization 
state in a beam of polarized light and said means for monitoring 
a polarization state in said polarized beam of light, there being 

„±_ — j— *i «...*.. 4- .-vwii -i^^,i^4-nUl a mr.nwc f av nnnf vaI 1 i nn ;a n 

pIKbeHL ctu i6asu One aujuDtauic m^c*n£> x.^a- w w*. v ± ^ 

ellipsometric phase angle between orthogonal components in a 
polarized beam of light, which adjustable means for controlling 
an ellipsometric phase angle, in use, allows sequentially setting 



a plurality of ell ipsometr ic phase angles between orthogonal 
components in a polarized beam of light which is caused by said 
ellipsometer system to interact with a material system, such that 
in use said ellipsometr ic phase angle can be set sequentially 
through a plurality of settings while ellipsometr ic data is 
obtained by said means for monitoring a polarization state in 
said polarized beam of light at at least two selected settings of 
said at least one adjustable means for controlling an 
ellipsometric phase angle; which obtained ell ipsometr ic data can 
be utilized in determination of PSI and DELTA values of an 
investigated material system, where said determination of said 
PSI and DELTA values includes compensating for the effects on 
said obtained ellipsometric data of said at least two selected 
setting of said at least one adjustable means for controlling an 
ellipsometric phase angle, on said obtained ellipsometric data; 

said ellipsometer system being further comprised of 
computational means which performs determination of investigated 
material system PSI and DELTA values, which computational means 
utilizes data obtained with said at least one adjustable means 
for controlling ellipsometric phase angle between orthogonal 
components in a polarized beam of light, being set to at least 
two selected settings, and which computational means performs 
compensation of the effects of said at least one adjustable means 
for controlling ellipsometric phase angle between orthogonal 
components, on said utilized ellipsometric data obtained at said 
at least two selected settings of said at least one adjustable 
means for controlling ellipsometric phase angle between 
orthogonal components, in determining material system PSI and 
DELTA values; 

b. placing a material system to be investigated into said 
ellipsometer system and causing a beam of polarized light from 
said means for setting at least one polarisation state in a beam 
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o£ polarized light to interact therewith and enter said means for 
monitoring a polarization state; 

c. adjusting said at least one adjustable means for controlling 
ellipsometric phase angle between said orthogonal components to 
be sequentially set to a plurality of settings while 
ellipsometric data is obtained by said means for monitoring a 
polarization state in said polarized beam of light at at least 
two selected settings from said plurality settings of said at 
least one adjustable means for controlling a value of 
ellipsometric phase angle between said orthogonal components; 

d. causing said computational means to determine investigated 
material system PSI and DELTA values by a method which performs 
compensation of the effects of said at least one adjustable means 
for controlling ellipsometric phase angle between orthogonal 
components in a polarized beam of light on said ellipsometric 
data obtained at said at least two selected settings of said at 
least one adjustable means for controlling ellipsometric phase 
angle between orthogonal components in a polarized beam of light 
which is caused to interact with a material system, in 
determining material system PSI and DELTA values; and 

e. optionally determining at least some of members of the 
group consisting of: 

the "Handedness"; 
Stokes Vector; and 

Jones and Mueller Matrix components; 

of said polarized beam of light and investigated material system. 

(It is noted that the computational means can be any computer 
system with sufficient memory and processing capability). 
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Said method can Include obtaining data comprising a 
plurality of relative magnitude ratios of orthogonal components 
and phase angles between orthogonal components are obtained, at 
least some o£ which plurality of ellipsometr ic relative magnitude 
ratios of orthogonal components and measured ellipsometr ic phase 
angles between orthogonal components correspond to sequential 
adjusted settings of ellipsometr ic relative magnitude ratios of 
orthogonal components present in said beam of polarized light, 
said sequential adjusted settings being effected by adjustment of 
at least one member of the group consisting of: 

said means for setting at least one polarisation state 

in a beam of polarized light; and 
said means for identifying a polarization state in said 

polarized beam of light; 

and in which said computational means is also caused to perform 
compensation of the effects of said sequential adjusted settings 
of ellipsometr ic relative magnitude ratios of orthogonal 
components present in said beam of polarized light which is 
caused by said ellipsometer system to interact with a material 
system, in determining investigated material system PSI and DELTA 
values. Said method can also involve obtaining a plurality of 
ellipsometr ic phase angles between orthogonal components which 
are effected at each sequential adjusted setting of ellipsometr ic 
relative magnitude ratio of orthogonal components present in said 
beam of polarized light which Is caused by said ellipsometer 
system to interact with a material system. Said 494 is disclosed 
as the at least one adjustable means for controlling 
ellipsometr ic phase angle between said orthogonal components to 
be sequentially set to a plurality of settings while 
ellipsometr ic data is obtained is a compensator which is 
typically utilized in practice of the present invention to 
determine "Handedness", that is, if a value for ( n ) which is 
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determined by that invention methodology is positive or negative* 

A recent Patent to Johs et al., No. 5,872,630, is also 
disclosed as it describes rotating compensator ellipsometer 
system and a method of evaluating material system 
characteristics, involving application of mathematical regression 
onto data acquired by application of electromagnetic beams 
thereto. Said 630 Patent Claims a spectroscopic rotating 
compensator system investigation system comprising a source of a 
polychromatic beam of electromagnetic radiation, a polarizer, a 
stage for supporting a material system, an analyzer, a dispersive 
optics and at least one detector system which contains a 
multiplicity of detector elements, said spectroscopic rotating 
compensator material system investigation system further 
comprising at least one compensator ( s ) positioned at a location 
selected from the group consisting of: 

before said stage for supporting a material system; and 

after said stage for supporting a material system; and 

both before and after said stage for supporting a 
material system; 

such that when said spectroscopic rotating compensator material 
system investigation system is used to investigate a material 
system present on said stage for supporting a material system, 
said analyzer and polarizer are maintained essentially fixed in 
position and at least one of said at least one compensator ( s ) is 
caused to continuously rotate while a polychromatic beam of 
electromagnetic radiation produced by said source of a 
polychromatic beam of electromagnetic radiation is caused to pass 
through said polarizer and said compensator (s ) , said 
polychromatic beam of electromagnetic radiation being also caused 
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to interact with said material system, pass through said analyzer 
and interact with said dispersive optics such that a multiplicity 
of essentially single wavelengths are caused to simultaneously 
enter a corresponding multiplicity of detector elements in said 
at least one detector system. Also Claimed is a regression 
based method of calibrating a spectroscopic rotating compensator 
material system investigation system comprising the steps of: 

a. providing a spectroscopic rotating compensator material 
system investigation system as described; 

b. developing a mathematical model of said spectroscopic 
rotating compensator material system investigation system which 
comprises as calibration parameter variables polarizer azimuthal 
angle orientation, present material system PSI, present material 
system DELTA, compensator azimuthal angle orientation(s ) , matrix 
components of said compensator (s } , and analyzer azimuthal angle 
orientation, which mathematical model is effectively a transfer 
function which enables calculation of electromagnetic beam 
intensity as a function of wavelength detected by a detector 
element, given intensity as a function of wavelength provided by 
said source of a polychromatic beam of electromagnetic radiation; 

c. causing a polychromatic beam of electromagnetic radiation 
produced by said source of a polychromatic beam of 
electromagnetic radiation, to pass through said polarizer, 
interact with a material system caused to be in the path thereof, 
pass through said analyzer, and interact with said dispersive 
optics such that a multiplicity of essentially single wavelengths 
are caused to simultaneously enter a corresponding multiplicity 
of detector elements in said at least one detector system, with 
said polychromatic beam of electromagnetic radiation also being 
caused to pass through said compensator ( s ) positioned at a 
location selected from the group consisting of: 
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before said stage for supporting a material system; 

after said stage for supporting a material system; and 

both before and after said stage for supporting a material 
system; 

d. obtaining an at least two dimensional data set of intensity 
values vs. wavelength and a parameter selected from the group 
consisting of: 

angle-of-incidence of said polychromatic beam of 
electromagnetic radiation with respect to a present 
material system; and 

azimuthal angle rotation of one element selected from 
the group consisting of: 

said polarizer; and 
said analyzer; 

over time, while at least one of said at least one compensator { s ) 
is caused to continuously rotate; 

e. performing a mathematical regression of said mathematical 
model onto said at least two dimensional data set, thereby 
evaluating calibration parameters in said mathematical model; 

said regression based calibration procedure evaluated calibration 
parameters serving to compensate said mathematical model for 
non-achromatic characteristics and non-idealities of said 
compensator (s), and for azimuthal angles of said polarizer. Said 
630 Patent is disclosed for general insight to ell ipsometr ic 
procedure . 
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Other Patents to Johs et al. are Nos. 5,666,201 and 
5,805,285. Said Patents, it is noted, describe, respectively, a 
spectroscopic ellipsometer and a spectroscopic spectrophotometer 
system, wherein a spectroscopic beam of electromagnetic radiation 
is caused to interact with a dispersive optics to the end that 
multiple orders of wavelengths are produced and detected. 

A Patent to Herzinger et al. No. 5,796,983 is identified as 
it describes the Kramers-Kroenig relationship as it applies to 
novel Dielectric Function Parametric model oscillator 
structures . 

A commonly owned co-pending Patent Application by Hilfiker 

and Herzinger, describes an ellipsometer/polar imeter , and/or 

spectrophotometer based method of determining values for 

anisotropic refractive indices (n ), (n ) and (n ) in, 

x y z 

respectively, "x", "y" and ,! z" directions in a material system, 

which ,f x", "y" and "z" directions have a determinable 

relationship to the orientation of an alignment surface of said 

material system. Said invention method, in its preferred 

embodiment, because of its sequence of steps and utilization of 

parameter values and/or data acquired in earlier steps in later 

steps, requires acquisition of relatively simple and easy to 

acquire data sets of only one (1) dimension where refractive 

indices (n ), (n ) and (n ) or differences therebetween (An ); 

x- y % " xy ' 

and (A n ) or (An ) are evaluated. However, variations of 
xz yz 

said preferred method allow for use of multiple dimension data 

sets in mathematical regression based determination of deviation 

from perpendicular to a material system alignment surface 

regarding the direction of an out-of-plane (n ) refractive index, 

z 

or where, for instance, separately identifiable steps are 
combined, or where it is desired to determine dispersion 
characteristics of one or more refractive indices which have 
previously been evaluated at a single wavelength. The essence 
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of the preferred present invention method can be recited as 
sequentially involving at least four basic steps, said four basic 
steps being: 



5 1*. 

Determining the precise orientation of a third orthogonal index 

(n ) of refraction which projects essentially perpendicular to an 
s 

alignment surface of a material system, said alignment surface 

being characterized by two essentially in-plane orthogonal, (ie. 

10 (n ) and (n )), indices of refraction; 
x y 

Causing a spectroscopic beam of electromagnetic radiation 
"I? originating from a source of a beam of electromagnetic radiation, 
15 ffii (which spectroscopic beam of electromagnetic radiation is 

preferably essentially linearly polarised and is comprised of a 
m plurality of wavelengths for which said material system is 
C essentially transparent), to approach said alignment surface of 
said material system along a locus which is essentially 
20 5 co-incident with the orientation of the third orthogonal index 

~£] (n ) of refraction, at least partially transmit through said 

IU z 

;£ J material system and enter a transmission detector to the end that 
a one-dimensional data set as a function of wavelength is 
acquired; and applying a computational means which is programmed 

2 5 

with a mathematical model for said material system to the end 

that a value for the difference between the in-plane (A n ) 

xy 

indices of refraction and the Euler angle orientations of said 
essentially in- 

30 are determined; 



essentially in-plane orthogonal (n^) (n^) indices of refraction 



3^ 

Causing a beam of electromagnetic radiation originating from said 
source of a beam of electromagnetic radiation, (which beam of 
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electromagnetic radiation is preferably essentially linearly 
polarized and is comprised of at least one wavelength) , to 
approach said alignment surface of said material system along a 
plurality of near normal angles-of -incidence to said material 
system alignment surface, interact with said material system and 
enter said transmission and/or reflection detector to the end 
that a one dimensional data set as a function of 
angle-of -incidence is acquired; and applying said computational 
means which is programmed with a mathematical model for said 
material system to the end that a value for the difference 
between at least one out~of~plane combination of indices of 
refraction selected from the group consisting of ^A n xz ^' 

and (An ) is determined; and 
^ yz 

Causing a beam of electromagnetic radiation originating from said 
source of a beam of electromagnetic radiation, (which beam of 
electromagnetic radiation is preferably essentially linearly 
polarised and is comprised of at least one wavelength) , to 
approach said alignment surface of said material system along a 
plurality of angles-of -incidence at near the Brewster condition 
to said material system alignment surface, interact with said 
material system and enter said reflection and/or transmission 
detector to the end that data as a function of angle-of- 
Incidence is acquired; and applying said computational means 
which is programmed with a mathematical model for said material 
system to the end that an absolute value for at least one index 
of refraction selected from the group consisting of ^ n x '/ ^ n y' 
and (n ) is/are directly determined, and, if desired or required, 

Z 

from previously determined values for (An ), (An ) and 

xy ^ xz 

(Aft ), determining absolute values for (n ), (n ) or (n ) not 
- yz J x y z 

directly evaluated* The present invention methodology is 
basically different in that it does not require precise 
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determination of the precise orientation of a third orthogonal 

index (n ) of refraction which projects essentially perpendicular 
z 

to an alignment surface of a material system. 

5 The present invention methodology is complementary to that in 

the pending Hilfiker et al. Application. 

Scientific articles which are of direct interest, and which 
are included herewithin by reference are: 

10 

"Optical Characterization Of Anisotropic Plastics", Hilfiker, 
Herzinger, Bungay, Woollam & Elman, Optical Interference 
Coatings, Op. Soc. Am, Tech. Dig. Series, Vol. 9, ((1938). 

1% "Characterization of Bi-Axially-Stretched Plastic Films By 
yy Generalized Ellipsometry" , Elman, Greener, Herzinger & Johs, Thin 
^ Solid Films, 313-314 (1998). 



"Phase And Microstructure Investigations Of Boron Nitride Thin 
2CE Films By Spectroscopic Ellipsometry In The Visible And Infrared 
;;{! Spectral Range", Franke, Schubert, Neumann, Tiwald, Thompson, 
|f| Woollam, Hahn a Richter, J. Appl. Phys . 82(6), (Sept. 1937). 

"Determination Of Optical Anisotropy In Calcite From Ultraviolet 
25 To Mid-Infrared By Generalized Ellipsometry", Thompson, DeVries, 
Tiwald, & Woollam, Thin Solid Films, 313-314 (1993). 

Additionally, an article by Johs, titled "Regression 
Calibration Method For Rotating Element Ellipsometers", which 
30 appeared in Thin Film Solids, Vol. 234 in 1993 is also identified 
as it predates the Chen et al . 350 Patent and describes an 
essentially similar approach to ellipsometer calibration. An 
article by Jellison Jr. titled "Data Analysis for Spectroscopic 
Ellipsometry", Thin Film Solids, 234, {13333 is identified as it 
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describes a method for determining the accuracy with which 
certain data points can be measured, which information allows 
adding a weighting factor to a curve fitting regression procedure 
as applied to a multiplicity of data points, said weighting 
factor serving to emphasize the effect of more accurate and 
precise data. A book by Azzam and Bashara titled "Ellipsometry 
and Polarized light" North-Holland, 1977 is disclosed and 
incorporated herein by reference for general theory. An article 
by Collins titled "Automated Rotating Element Ellipsometers : 
Calibration, Operation, and Real-Time Applications", Rev. Sci. 
Instrum. 61(8), August 1990 is identified as it provides insight 
into rotating element ellipsometers. An article by Kleim et al. 
titled "Systematic Errors in Rotat ing-Compensator Ellipsometry" 
published in J. Opt. Soc. Am. /Vol. 11, No. 9, Sept 1994 is 
identified as it describes calibration of rotating compensator 
ellipsometers. An Article by An and Collins titled "Waveform 
Analysis With Optical Multichannel Detectors: Applications for 
Rapid-Scan Spectroscopic Ell ipsometer " , Rev. Sci. Instrum., 62 
(3), August 1991 is also identified as it discusses effects such 
as Detection System Error Characterization, Stray Light, Image 
Persistence etc., and calibration thereof. Also disclosed are 
articles by Schubert et al . which describe "Generalized 
Ellipsometry". The first thereof is titled "Extension Of 
Rotating-Analyzer Ellipsometry To Generalized Ellipsometry: 
Determination Of The Dielectric Function Tensor From Uniaxial 
Ti02", J. Opt. Soc. Am. A. 13, (1996). This article describes 
methodology for evaluating all components of a Mueller Matrix 
dielectric function, but does not describe first diagonal izing 
said Mueller Matrix and evaluating components of only the three 
resulting complex dielectric functions or refractive indicies. 
The second such article is authored by M. Schubert alone and is 
titled "Polarization Dependent Parameters Of Arbitrary 
Anisotropic Homogeneous Epitaxial Systems", Phys . Rev. B 53, 
(1996). The third such article is titled "Generalized 
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Transmission Ellipsometry For Twisted Biaxial Dielectric Media: 
Application To Chiral Liquid Crystals", J. Opt. Soc. Am. A/Vol. 
13, No, 9 (1996). The Schubert et al . articles, it is noted, 
describe ellipsometer systems which have rotatable material 
5 system supporting stages, and primarily focus on investigation of 
material systems, utilizing said ellipsometer systems, which, (in 
general terms), material systems have isotropic "in-plane" 
refractive indices, (ie. where = n ). Further identified for 

authority regarding mathematical regression is a book titled 
10 Numerical Recipes in "C", 1988, Cambridge University Press. 

It is noted at this point that the present invention is 
O primarily a method of evaluating directions, and magnitudes of 
2{ Euler angles and real and imaginary components, of orthogonally 
i| related Kramers-Kroenig consistent complex refractive indicies in 
jfj optically thick material systems which present with an optical 
Ji axis oriented either in-plane or out-of -plane, with respect to an 
35 alignment surface thereof. The present invention methodology 
H works particularly well on optically thick material systems which 

ar e "uniaxial", in that corresponding real, and corresponding 
"X imaginary, components of at least two orthogonally related 
Q optically thick material system characterising diagonalized 
1 tensor : 
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30 complex dielectric functions or refractive indicies are of equal 
magnitude, said two entries being termed "ordinary", with the 
third entry being termed "extraordinary". (Note, the terminology 
"ordinary" and "extraordinary" have particular relevance, for 
instance, where an optically thick material system is composed 
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from a molecular lattice, the molecules of which have a 
longitudinal "c" axis and two orthogonally oriented axes an "a" 
plane. A dielectic function or refractive index aligned with the 
"c" axis, (which is generally the optical axis), is termed 
"extraordinary", while the dielectric functions or refractive 
indicies associated with the two orthogonally oriented axes can 
are termed "ordinary") ♦ 

Further, it is beneficial to note that randomly measured 
dielectric function or refractive index tensors which 
characterise a material system generally contain non-zero values 
in many of the six off-diagonal entries therein. Well known 
Matrix mathematical technique, however, provides that such 
matrices can be "diagonal ized" , thereby resulting in a Matrix 
which has non-zero entries only along a top-left to bottom-right 
diagonal. Matrix Diagonalization involves evaluating Eigenvalues 
of a given Matrix (M), then forming a Matrix (S) from columns of 
Matrices in which Eigenvalues are substituted thereinto, finding 
the inverse of the resulting Matrix (S" 1 ), and then multiplying 
(S)(M)(S" 1 ). The present invention works best where a 
Diagonalized Matrix in which the entries on the diagonal are two 
ordinary and one extraordinary complex refractive index or 
dielectric function, preferably wherein the two ordinary entries 
have equal magnitude real and equal magnitude imaginary 
components so as to be classified as representing a "uniaxial" 
optically thick material system. Where no diagonal entries have 
equal magnitude corresponding real, and equal magnitude 
corresponding imaginary components, a material system represented 
is termed "biaxial". Again, it is noted that the present 
invention methodology works particularly well for optically thick 
material systems which are properly classified as "uniaxial", 
where acquired data sets applied in dielectric function or 
refractive index complex component mathematical regression 
procedures are obtained utilising a single spectroscopic 
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electromagnetic beam angle-of -incidence, (of said spectroscopic 
electromagnetic beam to said alignment surface), and two 
optically thick material system rotation angle, (around a 
perpendicular to an alignment surface thereof), settings* Of 
course, data obtained at additional angles-of- incidence and 
rotational angles can also be utilized in a mathematical 
regression evaluation of dielectric function or complex 
refractive index real and imaginary component values, but such is 
not typically necessary where an optically thick material system 
is properly classified as "uniaxial", and even where off-diagnoal 
terms have non-zero values. 

Also, as further general background, it is to be appreciated 
that the terminology "Kramers-Kroenig refractive index 
consistency" refers to there being a mathematical integration 
relationship between real and imaginary components in a 
dielectric function or complex refractive index. That is, where 
one, (the real or imaginary), component is evaluated by a 
mathematical regression procedure, the other, (imaginary or real 
component, respectively), can be calculated therefrom by 
mathematical integration* In terms of dielectric constants this 
is demonstrated by the relationship: 



Where this relationship applies, (and it does in all practical 
cases), the required strength of a data set utilized in 
mathematical regression evaluation of real and imaginary 
components can be reduced as evaluating one complex component 
allows direct integration mediated calculation of the other, 
Alternatively, a sufficient data set to evaluate both real and 
imaginary components can be acquired and applied, and the results 
checked for Kramers-Kroenig consistency. {See Patent to 




o 
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Herzinger et al. No. 5,796,983 for additional insight to the 
Kramers-Kroenig relationship). As the immediately foregoing 
refers to both refractive indicies and dielectric constants, it 
is of benefit to point out the refractive indicies and dielectric 
constants are directly mathematically related by: 



given e = el + ie2 

then e 2 = el 2 - e2 2 +iele2 

and given n = SQRTCe 2 ) 

then n = SQRT{(el z - e2 2 ) + iele2) or n = n r + ik 
where n* = (el 2 - e2 2 ) 
k = iele2. 



where i = SQRT ( -1 ) . 



Thus, any reference to complex refractive index should be 
considered equivalent to reference to complex dielectric 
function, and vice-versa. 



It is also noted that a plane of incidence of a beam of 
electromagnetic radiation to an alignment surface of an optically 
thick material system includes its direction of propagation locus 
and a perpendicular to said alignment surface. Thus it can be 
stated that the angle-of -incidence of said beam of 
electromagnetic radiation, is the rotation angle between the 
locus of said beam of electromagnetic radiation and said 
alignment surface, taken around an axis defined as being both 
orthogonal to said direction of propagation locus and in the 
plane of the alignment surface of the optically thick material 
system. 



Further, it is noted that Euler angles describe the 
orientation of an optical axis of a material system with respect 
to the measurement geometry frame of reference . The Phi { ) 
Euler angle is that rotation angle around the normal to said 
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optically thick material system alignment surface necessary to 
align the optical axis projection in the alignment surface to the 
measurement frame of reference. The Theta {&) Euler angle is 
the angle of the optical axis in the resulting new coordinate 
system rotated around the resulting new "x"-axis direction; and 



the PSI ( Ui 5 Euler angle is the angle of the optical axis in this 



latest resulting new coordinate system rotated around its new 
"s"-axis . 

Finally, the terminology "optically thick" refers to a 
material system which provides a reflected beam of 
electromagnetic radiation caused to impinge thereupon, that is 
comprised primarily of components reflecting from the top surface 
thereof, and very little if any components reflected from a 
bottom surface thereof. Optically thick substrates can result 
from simply providing a very deep substrate in which reflections 
from a back-side are minimal because of attenuation, or can 
effectively be the result of making the substrate wedge shaped in 
cross-section or can effectively achieved by roughening the back 
side of the substrate. Optically thick dos not necesarily imply 
the material system is opaque at the wavelengths utilized. 

Even in view of the foregoing, a need remains for 
methodology for application in simultaneously evaluating Euler 
angles and direction, and magnitudes of complex components of 
orthogonally related dielectric function or refractive indicies 
in an optically thick material system, whether the optical axis 
thereof is oriented in-plane or out-of-plane with respect to an 
alignment surface thereof, and where said optically thick 
material system is characterised by three orthogonally related, 
Kramers-Kroenig consistent, refractive indicies. The method, 
while not limited to said application, is particularly applicable 
in the case of "uniaxial" optically thick material systems, 
wherein two of said three orthogonally related Kramers-Kroenig 
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consistent, optically thick material system characterizing 
diagonal ised tensor : 



£(E) = 



^ SC 

0 



0 0 

0 8„ 



complex dielectric function or refractive indicies have 
corresponding real components of equal magnitude and 
corresponding imaginary components of equal magnitude, (said 
entries being termed "ordinary", with the third entry being 
termed "extraordinary"). 
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DISCLOSURE OF THE INVENTION 



The preferred present invention method enables 
simultaneously evaluating mathematical model Euler Angles and 
5 parameters which describe directions, and magnitudes of real and 
imaginary components, of orthogonally related Kramers-Kroenig 
consistent complex dielectric functions or refractive indicies in 
optically thick material systems which present with an optical 
axis oriented in a selection from the group consisting of: 

10 

in-plane; and 
^ out -of -plane ; 

l_5g with respect to an alignment surface thereof, said optically 
yi thick material system being uniaxial in that corresponding real 
tfl and corresponding imaginary components of at least two 

orthogonally related optically thick material system 
0 characterizing diagonalised tensor; 

2(j|| 
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complex dielectric functions or refractive indicies are of equal 
magnitude, said method comprising, in any functional order, the 
steps of: 

30 a. providing an optically thick material system which presents 

with Kramers-Kroenig consistent complex dielectric functions or 
refractive indicies and with an optical axis oriented either 
in-plane or out-of-plane with respect to an alignment surface 



thereof, said optically thick material system being uniaxial in 
that corresponding real and corresponding imaginary components of 
at least two orthogonally related optically thick material system 
characterizing diagonalised tensor: 
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10 

complex dielectric functions or refractive indicies are of equal 
magnitude ; 

b. placing said optically thick material system into a system 
for directing at least one spectroscopic polarised 
electromagnetic beam(s) of radiation onto said alignment surface 
at at least one angle (s) of incidence removed from a normal to 
said alignment surface, said system for directing at least one 
spectroscopic polarised electromagnetic beam(s) of radiation onto 
said alignment surface at at least one angle(s) of incidence 
removed from a normal to said alignment surface comprising a 
reflection detector system; 

25 c. selecting at least one spectroscopic polarized 

electromagnetic beam(s) of radiation to at least partially 
comprise wavelengths for which said optically thick material 
system is non-transparent, and causing said at least one beam(s) 
of spectroscopic polarised electromagnetic radiation to impinge 

30 on said alignment surface of said optically thick material 

system, at at least one angle (s) of incidence removed from a 
normal to said alignment surface, in plane(s) of incidence which 
include the locus of said beam of spectroscopic polarised 
electromagnetic radiation and said normal to said alignment 
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surface, said at least one beam(s) of spectroscopic polarised 
electromagnetic radiation being caused to reflect from said 
alignment surface of said optically thick material system and 
into said reflection detector system; 

d. at said at least one angle(s) of incidence and at one 
rotation angle of said optically thick material system around 
said normal to the alignment surface thereof, obtaining 
reflection detector system mediated experimental reflection 
intensity data as a functions of wavelength and angle of 
incidence of said at least one beam of spectroscopic polarized 
electromagnetic radiation onto said optically thick material 
system alignment surface; 

e. from the data obtained in step d., determining Jones Matrix 
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components and if the off-diagonal components are essentially 
zero, rotate the optically thick material system about a 
perpendicular to the alignment surface thereof, and proceeding 
step f. and if not proceeding directly to step g.; 

£• at said at least one angle(s) of incidence and at one 
rotation angle of said optically thick material system around 
said normal to the alignment surface thereof, obtaining 
reflection detector system mediated experimental reflection 
intensity data as a functions of wavelength and angle of 
incidence of said at least one beam of spectroscopic polarized 
electromagnetic radiation onto said optically thick material 
system alignment surface; 
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g. simultaneously providing a mathematical model of said 
optically thick material system which includes as parameters 
therein at least: 

real and imaginary components for each of the 

orthogonally related, Kramers-Kroenig consistent 

tensor diagonal complex dielectric functions or refractive 

indicies; 

sufficient rotation about a normal to an alignment surface, 
and deviation from alignment with said alignment surface 
angle parameters to define the orientations of the 
orthogonally related, Kramers-Kroenig consistent tensor 
diagonal complex dielectric functions or refractive indicies, 
and orientation of the optical axis, with respect to the 
alignment surface; and 

Euler angles relating material system angles to a laboratory 
frame of reference; and 

h. via application of a mathematical regression technique, 
evaluating parameters in said mathematical model to provide a 
best-fit to acquired experimental data. 

A modified method of evaluating mathematical model 
parameters describing Euler angles and directions, and magnitudes 
of real and imaginary components, of orthogonally related 
Kramers-Kroenig consistent complex dielectric functions or 
refractive indicies in an optically thick material system which 
presents with an optical axis oriented in a selection from the 
group consisting of: 

in-plane; and 
out-of -plane ; 
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with respect to an alignment surface thereof. Said method 
basically comprises, in any functional order, the steps of: 

a. providing an optically thick material system which presents 
with orthogonally related Kramers-Kroenig consistent complex 
dielectric functions or refractive indicies and with an optical 
axis oriented either in-plane or out-of-plane with respect to an 
alignment surface thereof, said optically thick material system 
preferably, but not necessarily, being uniaxial in that 
corresponding real and corresponding imaginary components of at 
least two orthogonally related optically thick material system 
characterizing diagonalized tensor : 
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complex dielectric function or refractive indicies are of equal 
magnitude ; 

b. placing said optically thick material system into a system 
for directing at least one spectroscopic polarized 
electromagnetic beam(s) of radiation onto said alignment surface 
at at least one angle (s) of incidence removed from a normal to 
said alignment surface, said system for directing at least one 
spectroscopic polarized electromagnetic beam(s} of radiation onto 
said alignment surface at at least one angle (s) of incidence 
removed from a normal to said alignment surface comprising a 
reflection detector system; 

c* selecting at least one spectroscopic polarized 
electromagnetic beam(s) of radiation to at least partially 
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comprise wavelengths for which said optically thick material 
system is non-transparent, and causing said at least one beam(s) 
of spectroscopic polarized electromagnetic radiation to impinge 
on said alignment surface of said optically thick material 
system, at at least one angle(s) of incidence removed from a 
normal to said alignment surface, in plane(s) of incidence which 
include the locus of said beam of spectroscopic polarised 
electromagnetic radiation and said normal to said alignment 
surface, said at least one beam(s) of spectroscopic polarised 
electromagnetic radiation being caused to reflect from said 
alignment surface of said optically thick material system and 
into said reflection detector system; 

d. at said at least one angle (s) of incidence and at at least 
two rotation angles of said optically thick material system 
around said normal to the alignment surface thereof, obtaining 
reflection detector system mediated experimental reflection 
intensity data as a function of wavelength and angle of incidence 
of said at least one beam of spectroscopic polarised 
electromagnetic radiation onto said optically thick material 
system alignment surface; 

e. simultaneously providing a mathematical model of said 
optically thick material system which includes as parameters 
therein at least: 

real and imaginary components for each of the 
orthogonally related, Kramers-Kroenig consistent 
tensor diagonal complex dielectric function or refractive 
indicies; 

sufficient rotation about a normal to an alignment surface, 
and deviation from alignment with said alignment surface 
angle parameters to define the orientations of the 



orthogonally related, Kramers-Kroenig consistent tensor 
diagonal complex dielectric function or refractive indicies, 
and orientation of the optical axis, with respect to the 
alignment surface; and 

5 

Euler angles relating material system angles to a laboratory 
frame of reference; 

£• via application of a mathematical regression technique, 
10 evaluating parameters in said mathematical model in view of 

acquired experimental data. 

Pi Preferred present invention practice, where in-plane 

y3 projection directions of two orthogonally related complex 
15^ dielectric function or refractive indicies are known, involves 
C9 obtaining experimental reflection data at two rotation angles of 
^ said optically thick material system around said normal to the 
l| alignment surface thereof, said two rotation angles being 

appropriate to sequentially align a plane of incidence of said at 
20fK least one spectroscopic polarized electromagnetic beam of 

! f IJ radiation along the direction of an in-plane projection of two of 
if the orthogonally related complex dielectric function or 

refractive indicies which have equal magnitude corresponding real 
and equal magnitude corresponding imaginary components * 

25 

Alternatively, the present invention method can involve 
obtaining experimental reflection data at two rotation angles of 
said optically thick material system around said normal to the 
alignment surface thereof, where said two rotation angles are 
30 appropriate to align a plane of incidence of said at least one 

spectroscopic polarized electromagnetic beam of radiation along 
the direction of an in-plane projection of one of said two 
orthogonally related dielectric function or refractive indicies 
which have equal magnitude corresponding real and equal magnitude 
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corresponding imaginary components, but also involves obtaining 
experimental reflection data with said plane of incidence of said 
at least one beam of spectroscopic polarized electromagnetic 
radiation oriented along a direction rotationally between the 
directions of said in-plane projections of said two orthogonally 
related complex dielectric function or refractive indicies which 
have equal magnitude corresponding real and equal magnitude 
corresponding imaginary components. 

Further, particularly where an optical axis of an optically 
thick material system does not project in an alignment surface 
thereof, or perpendicular to said alignment surface, present 
invention practice can involve obtaining data sets at multiple, 
(eg. more than two), rotation angles of said optically thick 
material system around said normal to the alignment surface 
thereof, and said rotation angles can be selected to include, or 
not, orientations which place the plane of incidence of said at 
least one spectroscopic polarized electromagnetic beam of 
radiation along the in-plane projection direction of at least one 
of the two orthogonally related dielectric function or refractive 
indicies which have equal magnitude corresponding real and equal 
magnitude corresponding imaginary components. This approach is 
often necessary as in-plane projection directions of two 
orthogonally related complex dielectric function or refractive 
indicies are not known. 

It is also noted that in any version of the present 
invention methodology, the obtaining of experimental reflection 
intensity data as a function of wavelength and angle of incidence 
of said at least one beam of spectroscopic polarized 
electromagnetic radiation onto said optically thick material 
system alignment surface, can involve obtaining experimental data 
sets at two or more spectroscopic electromagnetic radiation beam 
angles of incidence and/or at two or more optically thick 
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material system rotation angles* 

It is noted that the mathematical regression technique 
preferably involves reduction of square-error. 

Present invention methodology can involve providing a 
transmission detector and obtaining transmission data at at least 
one angle(s) of incidence of said at least one beam of 
spectroscopic electromagnetic radiation removed from a normal to 
said alignment surface, This is not in contradiction to the 
material system "optically thick" designation as "optically 
thick" does not necessarily indicate opacity of a material system 
at a utilised wavelength. Rather it refers to a material system 
which does not significantly reflect incident electromagent ic 
radiation from a back-side thereof. 

Additionally, the present invention methodology can involve 
providing at least two optically thick material systems which 
each present with orthogonally related Kramer s-Kroenig consistent 
complex dielectric functions or refractive indicies and with an 
optical axis oriented either in-plane or out-of-plane with 
respect to an alignment surface thereof, each said optically 
thick material system being at least uniaxial in that 
corresponding real and corresponding imaginary components of at 
least two orthogonally related optically thick material system 
characterizing diagonalized tensor : 
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dielectric function or complex refractive indicies are of equal 
magnitude; said method further comprising, for each said 
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optically thick material system proceeding as described above, 
wherein said method further comprises, via application of a 
mathematical regression technique, simultaneously evaluating 
parameters in said mathematical models for each of the 
investigated optically thick material systems in view of acquired 
experimental data. 

It is also to be appreciated that the present invention 
method of simultaneously evaluating Euler angles and directions 
and magnitudes of real and imaginary components of orthogonally 
related Kramers-Kroenig consistent complex refractive indicies in 
an optically thick material system is particularly suited for 
application where the wavelengths selected at which said 
optically thick material system is non-transparent, are in the 
infrared range of 5 to 40 microns. 

Finally, it has been determined that a screening step is 
beneficial to practice of the present invention, said screening 
step being the evaluating of the components of a Jones Matrix, or 
equivalent elements of a Mueller Matrix. Only if said step 
provides Off-Diagonal terms which are not essentially zero (0.0), 
does one proceed. If said step returns Off -Diagonal terms which 
are essentially zero (0.0), the additional step of rotating the 
sample system being investigated is performed so that said 
Off-Diagonal terms become other than essentially zero (0.0). 

The present invention will be better understood by reference 
to the Detailed Description Section of this Specification, with 
reference to the Drawings. 
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SUMMARY OF THE INVENTION 



It is a primary purpose and/or objective of the present 
invention to provide a method of simultaneously evaluating 
mathematical model parameters describing Euler angles and 
directions, and magnitudes of orthogonally related complex 
components of dielectric function or refractive indicies in a 
optically thick material system which presents with an optical 
axis oriented either in-plane or out-of-plane with respect to an 
alignment surface thereof, said optically thick material system 
preferably, but not necessarily, being at least: 

uniaxial ; 

in that magnitudes of corresponding real, and magnitudes of 
corresponding imaginary, components of at least two orthogonally 
related optically thick material system characterizing 
diagonal i zed tensor : 



e(E) = 
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complex dielectric function or refractive indicies are equal. 

It is another purpose and/or objective of the present 
invention to, via acquisition and application of additional data, 
extend the present invention method of evaluating mathematical 
model parameters describing Euler angles and directions, and 
magnitudes of orthogonally related complex components of 
dielectric function or refractive indicies in an optically thick 
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material system which presents with an optical axis oriented 
either in-plane or out~of~plane with respect to an alignment 
surface thereof, to cases where said optically thick material 
system is: 

biaxial; 

in that no magnitudes of corresponding real, and magnitudes of 
corresponding imaginary, components of at least two orthogonally 
related optically thick material system characterizing 
diagonalized tensor: 
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complex dielectric function or refractive indicies are equal* 

It is another purpose and/or objective of the present 
invention to, via acquisition and application of additional data, 
extend the present invention method of evaluating mathematical 
model parameters describing Euler angles and directions, and 
magnitudes of orthogonally related complex components of 
dielectric function or refractive indicies in an optically thick 
material system which presents with an optical axis oriented 
either in-plane or out-of -plane with respect to an alignment 
surface thereof, to cases where said optically thick material 
system Is non-diagonal in that of f -diagonal components of a 
optically thick material system characterizing tensor are 
non-zero . 

It is another purpose and/or objective of the present 
invention to teach applying a preliminary screening step is 
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beneficial to practice of the present invention, said screening 
step being the evaluating of the components of a Jones Matrix, or 
equivalent elements of a Mueller Matrix, Only if said step 
provides Off-Diagonal terms which are not essentially zero (0.0), 
does one proceed. If said step returns Off-Diagonal terms which 
are essentially zero {0.05, the additional step of rotating the 
sample system being investigated is performed so that said 
Off -Diagonal terms become other than essentially zero (0-0}. 

Other purposes and/or objectives will become appearant by a 
reading of the Specification and Claims. 
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ERIEF DESCRIPTION OF THE DRAWINGS 



Fig. 1 shows a demonstrative Hexagonal Lattice structure with a 
"c" axis shown oriented vertically, and an "a" plane oriented 
orthogonal thereto * 

Fig. 2 shows that an optically thick material system arrangement 
realized by positioning a multiplicity of such hexagonal lattice 
structures horizontally lateral to one another. 

Fig* 3a demonstrates a side view of such an optically thick 
material system showing an extraordinary dielectric function 
projecting other than perpendicular to the surface (D), and one 
ordinary dielectric function projecting other than parallel to 
the surface (D) , 

Fig. 3b shows the Euler angles THETA (0), PHI ( (j) ) and 



Figs. 4a - 4d demonstrate Miller indicies of crystals cut (110), 
(111), (111) and (221) . 

Fig. 5a shows an optically thick material system and a beam of 
electromagnetic radiation oriented in a plane which directs it 
along the x-axis direction of a rectangular bulk T i° 2 * 

Fig. 5b shows the basic components of an el 1 ipsometer system. 

Figs. Sal 6c2 show PSI and DELTA data for (110) orientation 
Ti0 2 derived by regression onto data sets which were 
simultaneously applied in a wavelength by wavelength regression 
fit to provide the real and imaginary sc and pc dielectric 
function values shown in Figs. 6dl - 6e2. 
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Figs* 6dl and 6d2 show complex dielectric sc components (£ f ) and 
(6^) arrived at by simultaneous wavelength by wavelength 
regression onto data sets used to provide results shown in Figs. 
Gal - 6c2. 

Figs. 6el & 6e2. show complex dielectric pc components ( 6, ) and 
(6^) arrived at by simultaneous wavelength by wavelength 
regression onto data sets used to provide results shown in Figs. 
Sal - 6c2. 

Figs. 7al - 7c2 show PSI and DELTA data for (100) orientation 
TiO derived by regression onto data sets which were 

simultaneously applied in a wavelength by wavelength regression 
fit to provide the real and imaginary sc and pc dielectric 
function values shown in Figs. 7dl - 7e2. 

Figs* 7dl and 7d2 show complex dielectric sc components ( €r, ) and 
(6^) arrived at by simultaneous wavelength by wavelength 
regression onto data sets used to provide results shown in Figs. 
7al - 7c2. 

Figs. 7el & 7e2. show complex dielectric pc components ( ) and 

arrived at by simultaneous wavelength by wavelength 
regression onto data sets used to provide results shown in Figs. 
7al - 7c2. 

Figs. 8al - 8d2 show PSI and DELTA data for (111) orientation 
Ti0 2 derived by regression onto data sets which were 

simultaneously applied in a wavelength by wavelength regression 
fit to provide the real and imaginary sc and pc dielectric 
function values shown in Figs. 8el - 8f2. 

Figs. 8el and 8e2 show complex dielectric sc components (6,) and 
(£^3 arrived at by simultaneous wavelength by wavelength 
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regression onto data sets used to provide results shown in Figs. 
8al - 8d2. 

Figs* 8fl and 8£2 show complex dielectric sc components (£L f ) and 
(£ ) arrived at by simultaneous wavelength by wavelength 
regression onto data sets used to provide results shown in Figs, 
Sal - 8d2. 

Figs, 9a and 9b show real and imaginary pc dielectric constants 
for all three Ti0 2 optically thick material system orientations 

(110), (100) and (111) . 

Figs. 10a and 10b show real and imaginary sc dielectric constants 
for all three ^iO^ optically thick material system orientations 

(110), (100) and (111) . 
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DETAILED DESCRIPTION 

Turning now to Fig. 1, there is shown a demonstrative 
Hexagonal Lattice structure with a "c" axis shown oriented 
5 vertically, and an "a" plane oriented orthogonal thereto. Also 

shown is an extraordinary dielectric function or refractive index 
(£ e ) oriented along said "c" axis, and two orthogonal "ordinary" 
refractive indices (£ c , ) in a Surface (S). Fig. 2 shows that a 
simply arrived-at optically thick material system arrangement 
10 might be realized by positioning a multiplicity of such hexagonal 
lattice structures horizontally lateral to one another, and in 
f ^ such a simple optically thick material system the "c" axis, and 
yg coincident extraordinary dielectric function or refractive index 
2f are perpendicular to the surface (S), and the ordinary refractive 
ill indicies are "in-plane" with respect to the surface (S). In such 
a simply arrived-at optically thick material system the surface 
(S) can be considered an Alignment Surface, Proceeding, said 
§ simply arrived-at optically thick material system could be cut 
;? along an angle to said surface (S) provide an optically thick 
2DJ material system with an alignment surface which is not coincident 
"i! with the surface (S). That is, the lattice orientations in said 
jU optically thick material system would oriented so as to place 

said "c" axis and extraordinary refractive index direction other 
than perpendicular to a bulk physical optically thick material 
25 system alignment surface. It can be visualised, however, that 

were such a bulk optically thick material system as shown in Fig* 
2, cut along a diagonal (D) as shown in Fig. 1, that an bulk 
optically thick material system "alignment surface" would result 
in which the "c" axis and extraordinary refractive index 
30 direction would not project perpendicular thereto, and in which 
the ordinary refractive indicies would not project in-plane with 
said "alignment surface". Fig. 3a demonstrates a side view of 
such a resulting optically thick material system showing an 
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extraordinary dielectric function or refractive index projecting 
other than perpendicular to the surface (D), and one ordinary 
refractive index projecting other than parallel to the surface 
(D). Many such scenarios can be developed, and generally well 
known Miller indicies are utilized to designate how crystals are 
cut to arrive at various orientations. Figs. 4a - 4d demonstrate 
Miller indicies of crystals cut (110), (111), (111) and (221) 
respectively, and it can be determined that the approach to 
determining Miller indices involves identifying where a plane 
intercepts "x", "y" and "z" axes, then dividing that value into 
1. A (100) plane would then intercept the (, x" axis at 1.0 and be 
parallel to the "y"~"z" plane. 

For convenient reference, Fig. 3b is included to show the 



describe the orientation of an optical axis of a material system 
with respect to the measurement geometry frame of reference. The 
PHI { <p ) Euler angle is that rotation angle around the normal to 
said optically thick material system alignment surface necessary 
to align the optical axis projection in the alignment surface to 
the measurement frame of reference. The THETA (Q) Euler angle 
is the angle of the optical axis in the resulting new coordinate 
system rotated around the resulting new "x H -axls direction; and 
the PSI Euler angle is the angle of the optical axis in this 

latest resulting new coordinate system rotated around its new 
"z"-axis . 

In experimental ellipsometr ic work performed to demonstrate 
practice of the present invention methodology, rectangular shaped 
bulk TiO^ optically thick material systems oriented (100), (110) 

and (111), were investigated with electromagnetic radiation in 
the mid-IR range, (ie. wavelengths of between 5 and 40 microns)* 
For the (100) and (110) orientations the "c" axis, (which is the 
optical axis), is parallel to the surface probed by the 



Euler angles, THETA (0), PHI 




and PSI {tV). Euler angles 
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spectroscopic electromagnetic beam, and for the (111) 
orientation, the "c" axis is neither parallel or perpendicular to 
the surface probed by the spectroscopic electromagnetic beam. 

Fig. 5b shows the basic components of an ellipsometer 
system, A source of electromagnetic radiation (LS) is shown 
sequentially followed by a Polarizer (P), and optionally a 
Compensator (C). Also shown are a stage (STG) for supporting a 
Material System (MS), and both reflective (RM) and transmission 
(TM) mode Detector Systems (DET), which are preceded by an 
optional Compensator (C f ) (C ,f ), and a Polarizer (P). In use a 
beam (PPCLB) of polychromatic electromagnetic radiation is caused 
to exit said source of electromagnetic radiation (LS), have 
imposed upon it a state of polarization, interact with the 
optically thick material system (MS), then be detected. Change 
in polarization state in said beam (PPCLB) of polychromatic 
electromagnetic radiation in becoming (EPCLB), is determinative 
of Material System characteristics. 

Fig. 5a shows an optically thick material system and a beam 
of electromagnetic radiation oriented in a plane which directs it 
along the x-axis direction of a rectangular bulk TiO^ (110) 
and (100) orientations three orientations of said plane of 
incidence were utilized during data acquisition: 

1. "x-z" plane of incidence positioned along the rectangular 

Ti0 2 optically thick material system length (x) direction 

and perpendicular to the (y) direction; 

2. "x-z" plane of incidence positioned along the rectangular 

TiCU optically thick material system width (y) direction 

and perpendicular to the (x) direction; 

3. "x-z" plane of incidence positioned rotated between the 

rectangular TiO optically thick material system length 
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(x) and width (y) directions. 

When investigating the (110) orientation optically thick 
material system, when electromagnetic beam plane of incidence 
orientations #1 and #2 were utilized to place the plane of 
incidence along the optic axis and perpendicular thereto, no 
conversion between p and s polarizations occur. As a result 
measurement of the diagonal terms of a reflection Jones matrix 
using standard ell ipsometr ic measurement was sufficient. 
Respective ellipsometr ic PSI and DELTA Spectra, (where the angle 
of incidence of the electromagnetic beam was 72 degrees), for the 
orientation #1 are shown in Figs. 6al & 6a2, and for orientation 
#2 are shown in Figs. 6bl & 6b2. When orientation if 3 was 
utilized anisotropic ellipsometric PSI and DELTA data were 
acquired and are shown in Figs. 6cl and 6c2. Date acquired and 
demonstrated in Figs. 6al - 6c2 were simultaneously applied in a 
wavelength by wavelength regression fit to provide the real and 
imaginary sc and pc dielectric function values. Determined 
complex dielectric components ( ) and ( ) for sc are shown in 
Figs. 6dl & 6d2, and complex dielectric components ( ) and ( ) 
for pc are shown in Figs. Gel & 6e2. Said results were found to 
be Kramers-Kroenig consistent where values for "offset", Ec and 
A, (see Kramers-Kroenig equation in the Background Section), were 
evaluated to be: 

Best-fit parameters for Kramers-Kronig conversion of £i to Ei for (110) sample. 



Parameters 






for K-K fit 




offset. 


5.72±0.02 


5.47±0.02 


E c 


0.0226±0.0001 


0.0226±lxl0" 4 


A 


2.403±0.03 


3.83±0.03 



When investigating the (100) orientation optically thick 
material system, when orientation #3 was utilized anisotropic 
ellipsometric PSI and DELTA data were acquired and are shown in 
Figs. 7al and 7a2. As in the (110) case, in all #3, #1 and #2 
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orientations, the angle of incidence of the electromagnetic beam 
was 72 degrees* When electromagnetic beam plane of incidence 
orientations #1 and #2 were utilized to place the plane of 
incidence along the optic axis and perpendicular thereto, no 
conversion between p and s polarisations occur. As a result 
measurement of the diagonal terms of a reflection Jones matrix 
using standard ellipsometr ic measurement was sufficient. 
Respective ellipsometr ic PSI and DELTA Spectra, for the 
orientation #1 are shown in Figs. 7bl & 7b2, and for orientation 
It 2 are shown in Figs. 7cl & 7c2. The slight conversion between p 
and s was due to slight misalignment of the "c" axis in the (x-y) 
plane. Date acquired and demonstrated in Figs. 7al - 7c2 were 
simultaneously applied in a wavelength by wavelength regression 
fit to provide the real and imaginary sc and pc dielectric 
function values. Determined complex dielectric components {6r ( ) 
and (Gz) for sc are shown in Figs. 7dl & 7d2, and complex 
dielectric components (6 ( ) and ( e^) for pc are shown in Figs. 7el 
Sc 7e2. Said results were found to be Kramers-Kroenig consistent 
where values for "offset", Ec and A, (see Kramers-Kroenig 
equation in the Background Section), were evaluated to be: 
Best-fit parameters for Kramers-Kronig conversion of s 2 to si for (100) sample. 



Parameters for 


£l-sc 


£i-pc 


sample (100) K-K fit 




offset 


6.93±0.03 


5.89±0.01 


E c 


0.0234±0.0001 


0.0232±lxl0" 4 


A 


3.73±0.03 


2.4t±0.01 



In the case of the (111) sample, the "c" axis is not 
parallel to the (x-y) plane (ie. in-plane), or perpendicular 
thereto. Measurements were obtained with the plane o£ incidence 
oriented at four different locations, (eg. at -109, -18.8, -66.1 
and 21.54 degrees rotation of the optically thick material system 
around a perpendicular to an alignment surface thereof). 
Measurements at three of the orientations were obtained using 72 
degrees angle of incidence, and at one orientation, (ie. the 
21.54 degrees rotation), data was obtained using 30 and 55 
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degrees angles of incidence. Figs Sal, 8bl, 8cl, and 8dl, show 
corresponding ell ipsometr ic PSI results, and Figs. 8a2, 8d2, 8c2 
and 8d2 show corresponding ell ipsometr ic DELTA values. Date 
acquired and demonstrated in Figs. 8al - 8d2 were simultaneously 
applied in a wavelength by wavelength regression fit to provide 
the real and imaginary sc and pc dielectric function values. 
Determined complex dielectric components ( ) and ( ) for sc are 
shown in Figs. 8el & 3e2, and complex dielectric components ( ) 
and ( ) for pc are shown in Figs- 8fl & 8f2. Said results were 
found to be Kramers-Kr oenig consistent where values for "offset", 
Ec and A, {see Kramers-Kroenig equation in the Background 
Section), were evaluated to be: 

Best-fit parameters for Kramers-Kronig conversion of e 2 to 6] for (111) sample. 



1 Parameters for 
sample (11 1) K-K fit 


£l-sc 


£i-pc 


offset. 


5.53±0.04 


6.7±0.2 


E c 


0.0289±0.0001 


0.0244±0.0007 


A 


2.50±0.02 


3.2±0.2 



Figs. 9a and 9b show real and imaginary pc dielectric 
constants for all three TiC> 2 optically thick material system 

orientations (110), (100) and (111). Note that the fit is very 
good in all said cases, but the noise is present beyond 20 
microns for the (111) case 

Figs. 10a and 10b show real and imaginary sc dielectric 
constants for all three Ti0 2 optically thick material system 
orientations (110), (100) and (111). 

In all cases demonstrated, experimental data was fit using a 
wavelength by wavelength approach. 
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Having hereby disclosed the subject matter of the present 
invention, it should be obvious that many modifications, 
substitutions, and variations of the present invention are 
possible in view o£ the teachings. It is therefore to be 
understood that the invention may be practiced other than as 
specifically described, and should be limited in its breadth and 
scope only by the Claims. 
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